Direct observation of acoustic phonon mediated relaxation between coupled exciton 

states in a single quantum dot molecule 
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We probe acoustic phonon mediated relaxation between tunnel coupled exciton states in an in- 
dividual quantum dot molecule in which the inter-dot quantum coupling and energy separation 
between exciton states is continuously tuned using static electric field. Time resolved and tem- 
perature dependent optical spectroscopy are used to probe inter-level relaxation around the point 
of maximum coupling. The radiative lifetimes of the coupled excitonic states can be tuned from 
f-S ~ 2 ns to ~ 10 ns as the spatially direct and indirect character of the wavefunction is varied by 

I 1 detuning from resonance. Acoustic phonon mediated inter-level relaxation is shown to proceed over 

timescales comparable to the direct exciton radiative lifetime, indicative of a relaxation bottleneck 
for level spacings in the range AE ~ 3 — 6 meV. 



PACS numbers: 68.65.Hb, 78.55.Cr, 78.67.Hc 



The charge and spin degrees of freedom of localized charges in quantum dot (QD) nanostructures currently attract 
i-^J | significant interest due to their potential to form quantum bits for an inherently scalable quantum processor based 
£^ ■ on solid-state hardware. The implementation of such coherent devices using QDs can be traced to the strong three 
dimensional confinement that gives rise to relatively long coherence times for charge^ Q and spin0, Q excitations 
when compared to the time required for quantum state manipulation. Very recently, advances in QD fabrication 
have allowed the realization of self-assembled double dot nanostructures in which the quantum coupling between 
electronic states, and the spatial localization of the exciton wavefunctions, can be continuously varied by varying the 
potential applied to a metal gate contact. y| la, Q Such tunable QD-molcculcs (QDMs) may allow ultrafast coherent 
manipulation using electro-optical techniques and may exhibit conditional coherent dynamics. Q 
A main obstacle to the realization of coherent QDM devices is decoherence caused by coupling to phonons . TToL ITTl IT^ 
In pioneering work by Bockelmann|l3| phonon mediated relaxation between discrete electronic states in QDs was 
. predicted to occur over timescales comparable to the radiative lifetime, an effect commonly termed the phonon 
\Q ' bottleneck. However, experiments quickly revealed little or no evidence for bottleneck effects for level spacings of a few 
10s of meV. These findin gs w ere explained by the participation of multi phonon processes |14| . Coulomb scattering 
' |15| or polaron formation |l6fl In contrast, for level spacings of only a few meV, acoustic phonons have been shown to 
be the dominant source of decoherence for both single 0, and vertically coupled^, El El El E3 quantum dots. 
In this case the inter-level scattering rate can be strongly suppressed or enhanced, depending on the spatial extent 
of the carrier wavefunction. Therefore, experimental information on the phonon mediated coupling of the discrete 
states of a QDM is of fundamental importance for comparison with theory. 

In this paper we probe acoustic phonon mediated inter-level relaxation between molecular-like coupled single exciton 
(e + h) states in an individual self-assembled QDM. The inter-dot tunnel coupling is controlled by varying the axial 
electric field, enabling us to tune the coupled state energy splitting in the range AE ~ 3 - 6 meV and control the 
• i-H ! mixing between spatially direct (e + h localized in upper dot) and indirect (e in lower dot, h in upper dot) excitons. 
By combining continuous wave (CW), time resolved and temperature dependent spectroscopy we demonstrate the 
ability to control the exciton spontaneous emission lifetime and probe acoustic phonon mediated relaxation as the 
. . . i level spacing and spatial character of the exciton wavefunction is varied. 

The QDMs studied consisted of two layers of stacked InGaAs QDs with a nominal inter-dot separation of d = 10 
nm. The molecules were embedded in an n — i Schottky photodiode structure to control the static electric field (F) 
along the growth axis and Al shadow masks were defined to optically isolate single QDMs. A schematic of the device 
structure is presented in the inset of Fig. 1(a). Further details of the InGaAs-GaAs QDM epitaxial layer sequence 
and sample fabrication can be found elsewhere. 0, @ Time resolved studies were performed using a gated Si-APD 
and time-correlated single photon counting techniques that provided a temporal resolution of < 0.2 ns. Carriers were 
photogenerated by a mode-locked Ti:Sapphire laser tuned at 1.48 eV to excite the QDM wetting layer that delivers 
^2 ps duration pulses with an energy of 6 fJ. 
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FIG. 1: (color online) (a) PL spectra of a single QDM as a function of the electric field. (Inset) schematic of the n — i Schottky 
diodes studied. (b) Upper panel: Experimental (symbols) and calculated (dashed lines) exciton energies for Xa and Xb- Middle 
panel: Energy splitting between Xa and Xb and fitting (see text). Lower panel: Relative intensities Ia,b/{Ia + Ib) and the 
resonant field (dashed line). 



A series of time integrated photoluminescence (PL) spectra recorded as a function of electric field (F) are presented in 
the main panel of Fig. 1(a). A clear anti-crossing between two transitions, labeled Xa and X B , can clearly be observed. 
These two peaks, which will be the focus of this paper, arise due to the coupling of spatially direct and indirect neutral 
excitons as discussed in our previous publications. 0, Q Other peaks observed, even at the lowest excitation powers 
investi gate d (P ex ~ 0.1 W/cm 2 ) also arise from the same QDM and arise from negatively charged excitons as discussed 
in ref. fljj- At a low field (F < 14 kV/cm) Xb is predominantly an unmixed state with spatially direct character. 
Thus, it shifts weakly with field and has a high electron-hole overlap and, consequently, a large oscillator strength. In 
contrast, Xa is a state with predominantly indirect character over this low field range and, thus, shifts strongly with 
F and has a weak oscillator strength, jfj As the electric field increases, Xa and Xb shift closer to one another and 
become increasingly mixed until, for F cr n ~ 15.8 kV/cm, they anticross. At this point the electron component of the 
exciton wavefunction hybridizes over both dots to produce exciton states with bonding (Xb ) and anti-bonding (Xa ) 
character. The measured peak positions are compared with effective mass calculations of the coupled exciton states 
that include strain, quantum tunneling coupling, and Coulomb interaction for a model QDM with d = 11 nm. [(| 
The results of this calculation are shown by the dashed lines on Fig 1(b) (top panel), reproducing well both the field 
dependence of the transition energies of Xa and Xb and the coupling energy at F cr n marked by a vertical line. Below, 
we use these effective mass calculations to obtain the F-dependence of the relative oscillator strengths of Xa and 
Xb and compare them to the results of our time resolved and CW spectroscopy. A s expected for a coupled quantum 
system, the energy splitting between Xa and Xb is well represented by AE = (2AE e - e ) 2 + (S e - e • (F — F cr i t )) 2 
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FIG. 2: (color online) (a) Normalized PL decay curves for Xa and Xb, measured at F — F cr it = 2.6kV/cm. Dashed and the 
solid lines are fits with a mono-exponential function and a bi-exponential function, respectively, (b) The PL decay rates for 
Xa (D a ) (solid circles) and Xb (solid squares) obtained by the mono-exponential fittings. Total population decay rates for Xb 
(D B ) (open squares) based on the rate equation model discussed in the text are also shown. 



with (5 e _ e ■ (F — F cr it) being the detuning from the point of maximum coupling [Fig. 1(b) - middle panel]. 
The lower panel of Fig. 1(b) shows the relative intensities Ia,b/(Ia + Ib) of the upper and lower exciton branches 
as the resonance field is traversed. Clearly, the field at which the two exciton branches have equal relative intensities 
is larger than the critical electric field of F cr jt=15.8 kV/cm. This observation contrasts strongly with the simple 
expectation that the radiative rates and populations are similar for both states at -Fcrit ■ H3 Moreover, it indicates the 
presence of a relaxation process from Xa to Xb which is active over timescales comparable to the radiative lifetime. 
We note that the observed intensity ratio of I a /Ib ~ 1 : 5 does not reflect a thermal population distribution since 
only the lower exciton branch would be significantly populated [exp(— 2AF e _ e /fc,gT) < 10 -4 ]. 

In order to obtain further information about the nature of the relaxation process we performed time resolved 
spectroscopy on the two exciton branches as a function of electric field. Typical PL decay traces recorded for Xa and 
Xb are presented in Fig. 2(a) for an electric field detuned (F — F cr a — +2.6 kV/cm) from the critical field. Under 
these conditions Xb exhibits a much slower decay than Xa since the spatial overlap of the electron and the hole is 
much weaker for the mixed exciton state due to its stronger fraction of indirect character. 

Before entering into detailed analysis we present a summary of the decay lifetimes as a function of positive and 
negative detuning, either side of F crit . From fitting the observed decay transients using single exponential decay 
functions we obtain the decay rates for Xa and Xb presented in Fig. 2(b). With increasing the electric field the 
decay rate measured for Xa increases by a factor of ~ 3x whilst the Xb decay rate decreases by a similar factor. This 
distinct anti-correlation between the field dependent PL-decay rates recorded for Xa and Xb demonstrates clearly 
that the electric field tunes the direct-indirect character of the two mixed excitonic states. States with strong indirect 
character have lower decay rates, increasing as the direct admixture becomes larger. 

Surprisingly we find that Xa and Xb do not have the same decay rate at F cr u despite our simple expectation 
that the two states should carry similar oscillator strength. Instead, the observation of a different decay rate for Xa 
and Xb at F cr n strongly indicates the existence of a non-radiative relaxation channel from the antibonding (Xa) to 
the bonding state (Xb) with a rate comparable to the pure radiative lifetime. In order to quantify the relaxation 
rate and radiative lifetime, we consider a simplified rate equation model illustrated by the 3 level system presented 
in the inset of Fig. 3(a). In this model, the two neutral exciton levels Xa and Xb decay radiatively with electric 
field dependent rates R a and Rb, respectively. The non-radiative relaxation between the upper and lower exciton 
branches are taken into account by down (7) and up (7') scattering rates due to phonon emission. In addition, 
a phenomenological electron capture process with rate v is included due to the prominence of negatively charged 
excitons in the PL spectra presented in Fig. 1(a). The rate equations for the population of Xa (N a ) and Xb (Nj,) 
subject to the processes discussed above can be set in the form 
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FIG. 3: (color online) (a) Experimental (solid squares) and calculated (solid line) radiative recombination rates plotted as a 
function of F. (Inset) schematic processes for radiative recombinations of the exciton states Xa and Xb, relaxation between 
them, and electron capture, (b) Phonon relaxation rate 7 = 70 at 4.2 K from Xa to Xb obtained from the time-resolved PL 
and rate-equation model. 



^ = -R a N a - 1 N a + iN b - vN a + G, (1) 
= -R i N b -< J 'N b + <yN a -vN b + G, (2) 

where G is the neutral exciton generation rate which is the same for Xa and Xb- In order to simplify further our 
analysis, we note that 7' <C 7 = 70 at 4.2 K since the phonon population is negligible (< 10~ 4 ) and solve the rate 
equations to obtain the time dependence of N a and N b , viz. 

N a = exp(-D a t) (3) 
N b = (l + a)exp(-D b t)-aexp(-D a t). (4) 

Here D a = R a + v + 70 and D b = R b + v are the measured decay rates for Xa and X B and a = 7 /(D a — D b ). The 
decay rates extracted from a least squares fitting of mono- (N a ) and bi-exponential (N b ) transients are presented in 
Fig. 2(b) from which we estimate 70 ~ 0.6 ± 0.5 ns _1 . Here, the error is large as 70 is extracted from the relative 
amplitudes of the bi-exponential fit. Thus, we also estimated 70 by measuring the difference between D a and D b 
[Fig. 2(b)] at F cr it- This provides a direct measure of the inter-level relaxation rate since R a (F cr it) ~ Rb(F cr it) 
providing that the charged exciton formation rate is the same for both Xa and Xb ■ This approach yields 70 (-Fcrit) = 
D a (F cr it) — D b (F cr it) ~ 0.7 ± 0.3 ns _1 , in good agreement with the value obtained by fitting the bi-exponential decay 
transient. 

In order to investigate the F-dependence of 70 we need to estimate v, the electron capture rate. To do this we focus on 
the field range —1 kV/cm < (F—F crit ) < 3 kV/cm, since the charged excitons are pronounced at lower field and D b can 
be directly measured in this region. We expect that the recombination rate of the indirect exciton R b is much smaller 
than v in the high field range (F > 18 kV/cm), since Xb is much weaker than other lines arising from charged exciton 
recombination [Fig. 1(a)]. We then estimate v = 0.9 ns _1 from D b , constant in this narrow field range since the field 
dependence of the charged exciton intensity is much weaker than the neutral exciton transitions. Using this value for is, 
we deduce the radiative recombination rate R b . The results of this analysis is presented in Fig. 3(a) and compared with 
the calculated recombination rate obtained in the dipole approximation, i.e. R = n r e 2 E exc Ep \ (f e \fh)\ /2'RemQC'h 2 
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FIG. 4: (a) Experimental and (b) calculated relative intensities of Xa and Xb at 4.2 K and 45 K. 



where n r — 3.4 is the refractive index, E p = 25.7 eV is the Kane matrix element, E exc is the emission energy, and 
(f e \fh) is the overlap integral of the envelope functions obtained from the effective mass calculations used to fit the 
transition energies in Fig. 1(b). Fairly good quantitative agreement is obtained between these calculations and the 
values of Rb extracted from our rate equation analysis, strongly supporting the general validity of the analysis. The 
radiative lifetime can be tuned from ~ 2 to ~ 10 ns as the fraction of direct character in the mixed exciton state is 
tuned, in good quantitative agreement with previous reports for vertically stacked QDs.[2ll l22j From the extracted 
values of Rb and v we obtain the electric field dependence of 70 presented in Fig. 3(b). Typical values are of the order 
of 70 ~ 1 ns -1 for a level separation AE = 3 — 4 meV, consistent with recent theory for phonon mediated relaxation 
in QDMs with properties similar to those investigated here. [HE! 

At elevated temperatures both relaxation and activation processes determine the steady state populations N a and 
iVj, and, thus, the resulting intensities of Xa and Xb (I a and Is)- Fig- 4(a) shows the measured field dependence 
of the relative intensities at 4.2 K (AE > k B T, 7 > 7') and 45 K (AE ~ k B T, 7 ~ 7'). At 45 K, I A ~ I B for 
F ~ F cr i t as expected since R a ~ i?& at F cr u [Fig 3(b)] and the activation process results in N a ~ 2V&. This processes 
is reproduced by our rate equation model with single phonon scattering processes. Solving eqns. (1) and (2) in steady 
state we obtain the PL intensity ratio for a CW excitation at zero temperature, i.e. 



ia = i + {y + gyygb 

I B 1 + (v + 2j)/R a ' 



(5) 



At finite temperature, one-phonon emission (absorption) rates are given by 7(T) = 7o{A p ^(T) + 1} [^'(T) = 
JoNph(T)], where N p h(T) is the phonon occupation factor. The intensity ratios and their temperature dependence 
obtained by substituting the calculated R a , Rb and the estimated v into the above equation are presented in Fig. 
4(b), producing fairly good qualitative agreement with experiment. In particular, the shift of the equal intensity 
point towards F ~ F cr u is reproduced. The broad agreement between the intensity ratios obtained from the CW PL 
measurement and calculation strongly supports the validity of our rate equation analysis and the assumptions made. 
In summary, we have demonstrated the existence of inter-level phonon relaxation between the molecular-like exciton 
states in an individual self-assembled QDM. The relaxation rate is found to be comparable to the radiative recombi- 
nation rates of the direct excitons demonstrating the existence of a phonon bottleneck in the meV range. Furthermore, 
we tune the radiative recombination rates of a single molecule as the fraction of the direct exciton character is varied. 
The observation and the calculation of the thermal evolutions of the exciton populations suggests a dominance of 
one-phonon processes for the inter-level relaxation in the QDM, with an intrinsic rate ~ 0.5 — 1.5 ns -1 . 

This work was supported in part by the IT program (RR2002) from the MEXT (Japan) and by the DPG via 



6 



SFB631 (Germany). 



[1] A. Zrenner, E. Beham, S. Stufler, F. Findeis, M. Bichler, and G. Abstreiter, Nature (London) 418, 612 (2002). 
[2] T. Hayashi, T. Fujisawa, H. D. Cheong, Y. H. Jeong, and Y. Hirayama, Phys. Rev. Lett. 91, 226804 (2003). 
[3] M. Kroutvar, Y. Ducommun, D. Heiss, M. Bichler, D. Schuh, G. Abstreiter, and J. Finley, Nature (London) 432, 81 
(2004). 

[4] For a review, see, V. Cerletti, W. A. Coish, O. Gywat and D. Loss, Nanotechnology 16, R27 (2005). 

[5] H. J. Krenner, S. Stufler, M. Sabathil, E. C. Clark, P. Ester, M. Bichler, G. Abstreiter, J. J.Finley and A. Zrenner, New 
J. Phys. 7, 184 (2005). 

[6] H. J. Krenner, M. Sabathil, E. C. Clark, A. Kress, D. Schuh, M. Bichler, G. Abstreiter, and J. J. Finley, Phys. Rev. Lett. 
94, 057402 (2005). 

[7] G. Ortner, M. Bayer, Y. Lyanda-Geller, T. L. Reinecke, A. Kress, J. P. Reithmaier, and A. Forchel, Phys. Rev. Lett. 94, 
157401 (2005). 

[8] J. M. Villas-Boas, A. O. Govorov, and S. E. Ulloa, Phys. Rev. B 69, 125342, (2004). 
[9] J. I. Climente, A. Bertoni, G. Goldoni and E. Molinari, cond-mat/0604655 
[10] E. A. Muljarov, T. Takagahara, and R. Zimmermann, Phys. Rev Lett. 95, 177405 (2005). 

[11] P. Borri, W. Langbein, U. Woggon, M. Schwab, M. Bayer, S. Fafard, Z. Wasilewski, and P. Hawrylak, Phys. Rev. Lett. 
91, 267401 (2003). 

[12] P. Borri, W. Langbein, S. Schneider, U. Woggon, R. L. Sellin, D. Ouyang, and D. Bimberg Phys. Rev. Lett. 87, 157401 
(2001). 

[13] U. Bockelmann, G. Bastard, Phys. Rev. B 42, 8947 (1990). 
[14] T. Inoshita and H. Sakaki, Phys. Rev. B 46, 7260 (1992). 

[15] B. Ohnesorge, M. Albrecht, J. Oshinowo, A. Forchel and Y. Arakawa, Phys. Rev. B 54, 11532 (1996). 
[16] V. Preisler, T. Grange, R. Ferreira, L. A. de Vaulchier, Y. Guldner, F.J. Teran, M. Potemski, and A. Lemaitre, Phys. Rev. 
B 73, 075320 (2006). 

[17] T. Fujisawa, D. G. Austing, Y. Tokura, Y. Hirayama and S. Tarucha, Nature (London) 419, 278, (2002). 
[18] G. Ortner, R. Oulton, H. Kurtze, M. Schwab, D. R. Yakovlev, M. Bayer, S. Fafard, Z. Wasilewski, and P. Hawrylak Phys. 
Rev. B 72, 165353 (2005). 

[19] H. J. Krenner, E. C. Clark, T. Nakaoka, M. Bichler, C. Scheurer, G. Abstreiter, and J. J. Finley, submitted for publication, 
cond-mat/0604659 

[20] Our calculation gives the recombination- rate difference of 4 % between R a and Rb at F cr it- 
[21] C. Bardot, M. Schwab, M. Bayer, S. Fafard, Z. Wasilewski, and P. Hawrylak, Phys. Rev. B 72, 035314 (2005). 
[22] B. D. Gerardot, S. Strauf, M. J. A. de Dood, A. M. Bychkov, A. Badolato, K. Hennessy, E. L. Hu, D. Bouwmeester, and 
P. M. Petroff, Phys. Rev. Lett. 95, 137403 (2005). 



